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ABSTRACT 


A  technique  was  devised  for  obtaining  rate  laws  and  kinetic 
parameters  which  describe  the  thermal  degradation  of  plastics  from 
thermogravimetric  analysis  data.  The  method  is  based  on  the  inter- 
comparison  of  experiments  which  were  performed  at  different  linear 
rates  of  temperature  rise. 

The  method  was  applied  to  a  fiberglass  reinforced  phenol- 
formaldehyde  resin.  The  kinetic  equation  that  was  thus  derived 
was  checked  by  a  variety  of  techniques,  including  a  comparison 
with  constant  temperature  data  that  were  available  in  the  literature. 
The  tests  indicated  that  the  equation  was  applicable  to  a  major 
part  of  the  degradation  process. . 
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THE  KINETICS  OF  THERMAL  DEGRADATION  OF  CHARRING  PLASTICS* 

I.  GLASS  REINFORCED  PHENOL-FORMALDEHYDE 

Henry  L.  Friedman 

Introduction 

Kinetics  of  tnermal  degradation  of  solids  nave  been  evaluated 
from  tnermogravimetric  analysis  (TGA)  at  linear  rates  of  temperature 
rise  in  a  number  of  studies,  van  Krevelen,  van  Heerden,  and 
Huntjens^  studied  the  decomposition  of  coal.  Freeman  and  Carroll^ 
investigated  the  degradation  of  calcium  oxalate  monohydrate.  Their 
technique  was  applied  to  the  decomposition  of  a  styrenated  polyester 
by  Anderson  and  Freeman^,  and  was  modified  by  the  same  authors^  for 
polystyrene  and  polyethylene.  Doyle^  studied  the  degradation  of  a 
large  number  of  plastics.  Friedman^  investigated  the  decomposition 
of  a  glass  reinforced  phenol-formaldehyde.  A  more  sopnisticated 
approach  will  be  used  in  the  present  paper,  tnan  was  employed  for 
tne  earlier  work. 

In  all  of  the  earlier  studies,  calculations  were  based  on 
experiments  which  were  performed  at  a  single  rate  of  temperature  rise. 
In  the  present  paper,  kinetic  calculations  are  based  on  an  inter¬ 
comparison  of  experiments  which  were  carried  out  at  several  different 
rates  of  temperature  rise. 

*  This  work  was  performed  under  the  auspices  of  the  U.  S.  Air  Force 
Space  Systems  Division,  Contract  No.  AF  04(647) -269. 
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Experimental 


Fiberglass  reinforced  CTL  91-LD  Phenolic  (a  product  of 
Cincinnati  Testing  and  Research  Laboratory)  was  used  for  this 
study.  The  laminate  was  prepared  from  resin  impregnated  glass 
sheets  that  were  supplied  by  U,  S.  Polymeric  Chemicals,  Inc. 

The  sheets, which  originally  consisted  of  42  weight  per  cent  resin, 
were  B-staged  in  a  200°F  oven  for  one  hour*  They  were  then  cured 
in  a  225°F  press  at  100  p.s.i.  for  15  minutes,  followed  by  265°F 
at  100  p.s.i.  for  three  hours.  The  cured  laminate  was  allowed  to 
<  cool  to  room  temperature  in  the  press,  and  was  then  post  cured  in 
a  265°F  oven  for  16  hours. 

In  order  to  prepare  TGA  samples,  the  surface  layers  were  split 
off  and  discarded.  The  underlying  material  was  converted  to  thin 
flakes  by  end  milling  under  liquid  nitrogen.  The  flakes  were 
dried  in  a  desiccator  overnight,  and  were  further  dried  by  holding 
the  sample  in  a  porcelain  crucible  in  the  Chevenard  Thermobalance 
for  one  hour  at  100°C  just  before  TGA  was  allowed  to  commence. 
Samples  lost  less  than  one  per  cent  of  their  original  weight  during 
the  preconditioning  period. 

Four  200  mg  samples  of  the  flakes  were  subjected  to  TGA  at 
linear  heating  rates  of  50,  100,  180,  and  360°C/hr.  Degradation 
was  carried  out  in  a  stream  of  dry  nitrogen  to  maximum  temperature 
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of  900°C.  After  TGA  was  completed,  the  residue  from  each  test  was 
burned  in  air  at  1000°C  until  only  glass  remained.  Original  resin 
contents  were  found  to  range  from  35,0  to  36.1  weight  per  cent. 

Residual  weight  fractions  of  organic  material  were  calculated 
on  a  glass-free  basis,  and  are  shown  as  a  function  of  temperature 
in  Figure  1.  Note  that  the  curves  are  displaced  to  higher 
temperatures  with  increased  heating  rate,  as  would  be  predicted 
from  kinetic  theory.  Rates  of  weight  loss  are  plotted  as  a  function 
of  time  and  temperature  in  Figure  2.  Rates  were  measured  from 
enlarged  graphs  of  Figure  1  with  a  Gerber  Derivimeter.  All  of  the 
rate  curves  had  the  same  general  shape.  Slow  weight  loss  commenced 
almost  immediately  and  continued  to  above  300°C.  Between  300  and 
400°C,  the  rate  of  loss  of  weight  increased  as  the  temperatures 
were  elevated,  reached  a  maximum  in  the  vicinity  of  500°C  and  then 
decreased  for  the  remainder  of  the  heating  period.  Inspection  of 
Figure  2  shows  that  there  are  a  number  of  hills  and  valleys  in  the 
rate  curves.  Since  they  vary  markedly  from  curve  to  curve,  they  are 
thought  to  be  caused  by  scatter  of  the  TGA  data  rather  than  indications 
of  changes  in  the  degradation  process. 

The  early  weight  loss  accounted  for  less  than  10  per  cent  of 
the  total  loss,  and  appeared  to  be  associated  with  a  different 
process.  It  was,  therefore,  decided  to  eliminate  this  portion  from 
the  kinetic  analysis.  An  arbitrary  initial  weight  of  plastic  was 
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chosen  from  the  TGA  data  at  a  temperature  between  300  and  400°C 
where  examination  of  the  rate  curves  showed  that  the  major  weight 
loss  process  appeared  to  begin.  Although  it  was  likely  that  the 
main  degradation  process  started  at  somewhat  lower  temperatures, 
and  that  the  early  degassing  process  continued  to  higher  temperatures, 
it  was  felt  that  the  errors  associated  with  rate  measurements  would 
be  far  more  serious  than  errors  caused  by  a  reasonable  estimate  of 
initial  weight.  Figure  3  shows  residual  weight  fractions  of 
organic  material,  with  the  early  weight  loss  process  eliminated. 
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Kinetic  Analysis 


The  following  very  general  kinetic  equation  was  assumed  to 
hold  for  this  material. 

-  i-  g  -  A  e'AE/RT  £  (v/Wo)  (1) 

where  w  ■  weight  of  organic  material. 
wQ  -  original  weight  of  plastic, 
t  -  time,  hrs. 

A  ■  pre-exponential  factor  of  rate  constant,  hr.”^. 

/\E  «  activation  energy  of  rate  constant,  cal./mole. 

R  »  gas  constant,  1.987  cal. /°K  mole. 

T  -  absolute  temperature,  °K. 
f  (w/w0)  ■  a  function  of  the  weight  of  organic  material. 

Taking  logarithms  of  (1)  gives 

ln  (-  °  In  A  +  In  f  (w/w0)  -  (2) 

f  (w/wQ)  is  assumed  to  be  constant  for  constant  values  of  w/wQ.  This 
is  comparable  to  assuming  that  the  chemistry  of  the  process  is 
independent  of  temperature,  and  is  dependent  only  on  the  instantaneous 
weight  of  the  organic  material.  A  dozen  values  of  w/wQ  were  selected, 
ranging  from  0.675  to  0,950,  at  equal  intervals.  Values  of  (l/wQ) 
(dw/dt)  and  T  were  determined  for  each  w/w0  for  each  TGA  experiment. 
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Plots  of  In  (-l/wQ)  (dw/dt)  vs.  1/T  are  shown  for  each  value  of 
w/w0  in  Figure  4.  The  slope  of  each  line  is  equal  to  -AE/R,  while 
the  intercept  is  In  J_  A  x  f  (w/wQ)_7.  The  lines  of  Figure  4  were 
obtained  from  least  squares  treatments  of  the  data.  Z\E's  and 
In  /  A  x  f  (w/w0)_/'s  are  plotted  as  functions  of  w/wQ  in  Figure  5. 
Z\E  varies  from  50  to  72  kcal/mole.  The  range  of  variation  that 
is  associated  with  each  point  represents  only  the  probable  error 
of  the  points  based  on  the  least  squares  treatment.  The  real 
errors  are  larger.  Note,  however,  that  trends  in  the  values  of 
Z\E  and  In  /  A  x  f  (w/wQ)_/  are  identical.  This  indicates  that 
a  positive  error  in  A,  which  would  make  the  reaction  appear  to  go 
faster,  is  compensated  for  by  a  positive  error  in  Ae,  which  would 
tend  to  make  the  reaction  appear  to  go  slower.  Negative  errors 
would  balance  in  the  reverse  way.  Thus  it  appears  that  the  kinetics 
of  the  overall  process  are  consistent. 

The  12  different  AE 1  s  were  averaged  to  give  57.7  kcal.  Experi¬ 
mental  parameters  and  the  Ae  were  substituted  in  (2)  to  give  48 
values  of  In  j_  A  x  f  (w/w0)_/.  Average  In  /  A  x  f  (w/w0)_/'s  are 
plotted  as  a  function  of  w/wQ  in  Figure  6.  The  range  associated 
with  each  point  is  the  standard  deviation  of  the  average  of  four 
values.  The  smooth  curve  in  Figure  6  tended  to  confirm  that  the 
kinetics  of  the  overall  process  were  consistent,  f  (w/wQ)  was 
assumed  to  have  the  form 
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f  (w/w0)  -  (^-^)n  (3) 

wo 

where  Wf  «  final  weight  of  char. 

n  -  kinetic  order  of  the  reaction, 
w  -  wf  is  a  measure  of  the  amount  of  material  which  is  available 
for  decomposition  at  any  instant  of  time,  and  is  a  convenient 
mathematical  form  which  allows  the  degradation  to  be  completed  when 
all  of  the  original  plastic  is  converted  to  char.  Equation  (3),  or 
very  similar  forms,  can  be  derived  from  theory  for  processes  where 
a  solid  reactant  is  converted  to  a  constant  ratio  of  gaseous  and 
solid  products.  Multiplying  both  sides  of  (3)  by  A  and  taking 
logarithms , 

In  J_  A  x  f  (w/wQ)_/  -  In  A  +  n  In  (- - — )  (4) 

wo 

If  the  assumed  form  of  f  (w/wQ)  is  correct,  then  a  plot  of 
In  /  A  x  f  (w/w0)_/  vs .  In  [_  (w  -  wf)/w0_/  should  give  a  straight 
line  where  the  slope  is  n  and  the  intercept  is  In  A.  This  type 
of  plot  is  shown  in  Figure  7.  A  straight  line  resulted  for  all 
values  of  the  abscissa  below  -  1.3  (equivalent  to  w/w0  below  0.9). 
The  positive  deviation  at  higher  weights  could  be  attributed  to 
contributions  of  the  early  weight  loss  process.  In  calculating 
values  of  (w  -  wf)/wQ,  it  was  necessary  to  assume  a  value  of  Wf. 
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Examination  of  Figure  3  indicated  that  decomposition  was  essentially 
completed  at  900°C,  so  the  convenient  average  value  of  W£/wQ  ■  0.61 
was  selected.  Calculations  based  on  Figure  7  showed  that  n  ■  5.0 
and  A  -  5.98  x  10^®  hr.”^.  The  form  of  the  overall  kinetic  equation 
is 

_  1_  dw  _  i  w  -  wf. 5  (5) 

w0  dt  k  <-*— > 

where  k  -  A  e“  ^E/RT f  hr."^-. 
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Tests  of  the  Kinetic  Equation 
Three  different  tests  were  made  on  the  applicability  of 
equation  (5).  The  first  was  an  examination  of  Arrhenius  plots  for 
the  individual  runs,  k's  were  calculated  from  (5)  by  using  the 
values  of  the  otner  terms  which  were  calculated  earlier  in  the 
paper.  Wf/wQ  was  assumed  to  be  0.61  for  all  but  the  100°C/hr. 
run,  where  the  value  of  0.59  was  employed.  Arrhenius  plots  of  the 
four  runs  are  shown  in  Figure  8.  The  lines  are  based  on  visual 
inspection  of  the  data.  The  </\E's  and  A's  which  were  derived  from 
the  lines  are  listed  in  Table  1.  Note  that  the  values  of  k  are 
higher  than  the  Arrhenius  lines  at  the  lowest  temperatures.  This 
is  thought  to  be  due  to  the  early  weight  loss  process,  and  is 
comparable  to  the  deviations  from  linearity  in  Figure  6.  The 
method  of  least  squares  was  applied  to  all  the  data  of  Figure  8, 
except  for  the  two  lowest  temperatures  of  each  TGA  run.  /\E  was 

18  ~  i 

found  to  be  55  kcal,  while  A  is  10  hr.  .  These  values  will  be 

used  for  all  further  calculations. 

The  second  test  is  based  on  the  fact  that  at  the  maximum 

2  2 

(-dw/dt),  the  second  derivative  (-d  w/dt  )  is  equal  to  zero. 

Since  the  heating  rate  is  linear 


dT 

dt 


where =  heating  rate,  °C/hr. 
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Substituting  (6)  in  (5)  and  rearranging  gives 


-  <*-  *  ,  e-AE/RT  ,  .  y 

It  Tv?  (  f) 


Taking  derivatives  of  (7) 

.  d^w  Aw-  wf  4  -AE/RT 
d?  (~ - } 


5  dw  (  .  )  £& 

dT  v  rx2 


For  (-d^/dT2)  -  0  at  the  maximum  rate 

5  ~  +  (w  -  wf)  ^2  "  0 
dT  1  RTn? 

where  Tm  ■  temperature  at  maximum  rate,  °K. 
(9)  may  be  arranged  to 

AE  -  5RTm2  (-dw/dT) 

W  -  Wf 


(7) 


(8) 


(9) 


(10) 


The  four  values  of  /\E  which  were  calculated  from  (10)  are  listed  in 
Table  1. 

The  third  test  was  to  derive  the  individual  TGA  curves  from 
equation  (5).  (5)  was  rearranged  for  integration  to  give 


where  wQ 


dw 

(w  -  w0)3 


7^? 


e-AE/RT  dT 


(ID 


wo  'Ao 

■  the  arbitrary  choice  of  initial  weight  of  plastic. 
-  temperature  at  wQ,  °K. 
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(12) 


Integration  and  solution  for  w  gives 


w  -  wf  + 


where  v  -  AE/RT  and  v0  -  /SE/RT0. 

Solutions  for  the  integrals  which  appear  in  (12)  are  tabulated  in 
7  8  8 

several  sources  *  .  Reference  was  used  for  this  work.  Figure  9 
shows  a  typical  thermogram  that  was  derived  by  this  technique. 

Note  that  in  the  lower  temperature  range  the  experimental  curve 
shows  a  higher  weight  loss  than  the  theoretical  curve.  This  was 
expected,  since  the  kinetic  analysis  and  Arrhenius  plots  had  both 
shown  that  the  early  weight  loss  process  continues  after  temperatures 
which  are  greater  than  T0. 

A  plot  of  actual  vs.  calculated  w/w0's  for  all  TGA  curves  is 
shown  in  Figure  10.  Note  that  the  deviations  from  the  45°  line  at 
the  higher  mass  end  are  rather  systematic.  The  curve  that  is  shown 
may  be  used  to  correct  calculated  values  of  w  in  the  region  where 
the  early  weight  loss  process  is  important. 
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Discussion 


The  various  tests  of  equation  (5)  showed  that  the  kinetic 
analysis  of  CTL  91-LD  phenolic  degradation  that  was  formulated  in 
this  paper  is  reasonable  for  a  large  part  of  the  decomposition 
process.  The  main  deviation  between  experiment  and  analysis 
occurred  during  early  degradation,  where  the  actual  weight  loss 
exceeded  the  calculated  loss.  The  earliest  weight  loss  can  probably 
be  attributed  to  desorption  of  adsorbed  gases.  The  resin  was 
probably  not  completely  cross-linked  at  the  point  where  TGA  started. 
As  heating  continued,  further  cross-linking  probably  occurred,  and 
further  weight  loss  cou^d  be  attributed  to  the  evolution  of  some 
condensation  products.  Since  degradation  probably  started  before 
cross-linking  was  completed,  there  is  probably  no  simple  value 
for  w0,  in  the  way  that  it  was  used.  Examination  of  Figure  9  shows 
that  a  slightly  lower  wQ  would  make  the  theoretical  curve  fit  the 
experiment  better.  It  is  altogether  reasonable,  when  one  considers 
the  argument  that  was  presented  above,  that  the  effective  wQ  be 
lower  than  the  chosen  value.  It  is  also  not  surprising  that  the 
theoretical  weight  at  900°C  is  slightly  higher  than  the  experimental 
figure.  Adjustments  could  probably  be  made  in  the  choice  of  both 
wQ  and  Wf,  so  that  very  accurate  theoretical  TGA  curves  would  result, 
except  for  the  very  earliest  weight  loss  processes. 
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A  fifth  order  weight  dependence  is  unusual  for  polymer 

degradation.  Of  course,  order  does  not  mean  the  same  thing  in 

this  type  of  process  as  it  does  for  a  gas  or  a  solution.  Since 

this  resin  is  a  highly  cross-linked  material,  it  is  likely  that 

the  high  order  was  associated  with  a  rather  complex  mechanism 

which  occurred  in  the  solid  state.  Note,  however,  that  a  pre- 

14  -1 

exponential  factor  of  about  3  x  10  sec.  and  an  activation 
energy  of  55  kcal.  are  not  unusual  for  some  thermoplastics . 

Inspection  of  the  low  temperature  end  of  the  TGA  curves 
showed  that  very  low  activation  energies  could  be  attributed  to 
that  part  of  the  decomposition.  In  the  earlier  kinetic  analysis^, 
this  author  studied  TGA  experiments  which  were  performed  only  at 
180°C/hr.  on  a  sample  where  the  resin  content  had  not  been 
determined.  Since  that  analysis  covered  temperatures  ranging 
between  404  and  524°C,  it  is  reasonable  that  an  activation  energy 
of  15  kcal.  was  calculated.  Madorsky  and  Straus^ investigated 
the  thermal  degradation  of  this  material  in  vacuum  and  at  constant 
temperature.  All  of  their  rate  measurements  were  concerned  with 
early  degradation  processes.  They  estimated  an  activation  energy 
of  18  kcal.  from  their  data.  It  is  felt  that  none  of  the  work 
concerned  with  early  degradation  is  definitive,  and  that  this  region 
requires  further  study. 
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Madorsky  and  Straus  carried  out  a  few  experiments  to 
comparatively  high  extents  of  reaction,  and  it  is  possible  to 
compare  their  data  with  that  of  the  present  work.  Total  percent 
of  volatilization  was  reported  as  a  consequence  of  heating  for  a 
given  length  of  time  at  420,  500,  and  800°C.  To  compare  data, 
equation  (5)  was  integrated  to  give 

w0  (w0  -  Wf) 

y  b  -}■  *  "  *  . * — . — 1  . . .  .—  -I  -  -  i  .  .  (13) 

/■»„*  +  4  X  1018  t  e-55,000/RI  (wo  .  „f)4_7°-25 

These  workers  pyrolyzed  a  sample  of  phenolic  at  1200°C^8.  Elemental 
analysis  showed  that  99.2%  of  the  residue  was  carbon  when  the  residue 
weight  was  51.6%  of  that  of  the  original  sample.  wQ  was  assumed  to 
be  unity  ani  Wf  was  assumed  to  be  0.516.  Although  equation  (13) 
ignored  early  weight  loss  and  w0  was  defined  on  a  different  basis 
than  that  by  which  the  equation  was  derived,  it  was  felt  that  the 
calculation  would  still  be  meaningful.  Wf  is  probably  the  most 
important  parameter  for  this  calculation  so  long  as  the  extent  of 
reaction  is  great  enough  for  most  or  all  of  the  early  weight  loss 
gases  to  be  evolved.  The  three  runs  that  were  chosen  seemed  to 
qualify  reasonably  well  on  this  basis.  The  results  of  comparison 
are  listed  in  Table  2.  The  agreement  between  the  results  is 
excellent  and  lends  further  support  to  the  present  kinetic  analysis. 

The  kinetic  analysis  technique  that  was  employed  in  this 
paper  has  several  advantages.  By  its  use  it  was  possible  to 
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calculate  an  activation  energy  for  the  main  degradation  process 
without  any  knowledge  of  the  form  of  the  kinecxc  equation.  It 
was  also  useful  for  the  purpose  of  deriving  a  semi-empirical  rate 
law,  which  would  be  difficult  to  arrive  at  by  other  methods.  This 
method  could  also  serve  as  a  check  for  systems  where  the  rate  laws 
seem  to  be  so  simple  that  TGA  experiments  would  only  be  carried 
out  at  one  rate  of  heating  under  normal  circumstances. 
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A  technique  was  devised  for  obtaining  rate  laws  and  kinetic 
parameters  which  describe  the  thermal  degradation  of  plastics 
from  TGA  data.  The  method  is  based  on  the  intercomparison  of 
experiments  which  were  performed  at  different  linear  rates  of 
heating.  By  this  method  it  is  possible  to  determine  the  activation 
energy  of  certain  processes  without  knowing  the  form  of  the  kinetic 
equation. 

This  technique  was  applied  to  fiberglass  reinforced  CTL  91-LD 
Phenolic,  where  the  rate  law 

-  i_  dw  =  1018  e-55»000/RT  -  wfx5,  hr."1 

w0  dt  1  w0  > 

was  found  to  apply  to  a  major  part  of  the  degradation.  The  equation 
was  successfully  tested  by  several  techniques,  including  a  comparison 
with  constant  temperature  data  that  were  available  in  the  literature. 
The  activation  energy  was  thought  to  be  correct  to  within  10  kcal. 
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*Arrhenius  plot  average  values  were  calculated  from  examination  of  all  data  in  Figure 


TABLE  2 


Comparison  Between  Results  of  the  Present  Work  with 
Experiments  of  Madorsky  and  Straus 


Temperature, 

°C 

Reference 

Heating 

Time, 

Minutes 

w/w0 

observed 

w/w0 

calculated 

420 

10 

130 

0.838 

0.876 

500 

9 

30 

0.72 

0.717 

800 

9 

5 

0.53 

0.542 
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Figure  2.  Derivatives  of  TGA  Curves 


Figure  3.  TGA 
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Figure  4.  Arrhenius  Type  PL 
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Figure  7.  Plot  for  the  Determination  of  A  and  Kinetic  Order 


W/Wo-  ACTUAL 


1.0 


tween  Experiment  and  Theory. 


Figure  10.  Comparison  Be 


GENERAL®  ELECTRIC 


SPACE  SCIENCES  LABORATORY 

MISSILE  AND  SPACE  VEHICLE  DEPARTMENT 


TECHNICAL  INFORMATION  SERIES 


SUBJECT  CLASSIFICATION 


Henry  L. 
Friedman 


Plastic  Degradation 
Kinetics 


R61SD145 
August  1961 


TITlE  THE  KINETICS  OF  THERMAL  DEGRADATION  OF 
CHARRING  PLASTICS.  I.  GLASS  REINFORCED  PHENOL- 
FORHALDEHYDE _ _  _  _ ______ 

ABSTRACT 

A  tecnnique  was  devised  for  obtaining  rate 
laws  and  kinetic  parameters  which  describe  the 
tnermal  degradation  of  plastics  from  thermo- 
gravimetric  analysis  data.  The  method  is  based 
on  the  inter-comparison  of  experiments  which 


GOV.  CLASS 

Unclassified 


REPRODUCIBLE  COPY  FILED  AT 
G.  E.  TECHNICAL  INFORMATION  CENTFR 
3198  CHESTNUT  STREET 
PHILADELPHIA,  PENNA. 


Abstract  (Continued) 

were  performed  at  different  linear  rates  of 
temperature  rise. 

The  method  was  applied  to  a  fiberglass 
reinforced  phenol-formaldehyde  resin.  The 
kinetic  equation  that  was  tnus  derived  was 
checked  by  a  variety  of  techniques,  including  a 
comparison  with  constant  temperature  data  that 
were  available  in  the  literature.  The  tests 
indicated  that  the  equation  was  applicable  to  a 
major  part  of  the  degradation  process. 
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